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NOTATION 


A* 

a 

b 


-N 


^Pc 


J yo 


ya+ 


ya- 


zo 


za 


■ a+ 


°Pa- 


1 , 2,3 


relative area of a sector of an elliptical quadrant (sketch (a)) 
semiminor axis of an ellipse 
semimajor axis of an ellipse 

chord-force coefficient based on base area of elliptic cone 
drag-force coefficient based on base area of elliptic cone 
lift-force coefficient based on base area of elliptic cone 
normal-force coefficient based on base area of elliptic cone 
surface pressure coefficient 

surface pressure coefficient on circular cone at zero angle of 
attack 

surface pressure coefficient on two-dimensional wedge 

surface pressure coefficient on elliptic cone at zero angle of 
attack in vertical plane of symmetry 

windward surface pressure coefficient on elliptic cone at angle of 
attack in vertical plane of symmetry 

leeward surface pressure coefficient on elliptic cone at angle of 
attack in vertical plane of symmetry 

surface pressure coefficient on elliptic cone at zero angle of 
attack in horizontal plane of symmetry 

surface pressure coefficient on elliptic cone at angle of attack in 
horizontal plane of symmetry 

windward surface pressure coefficient on elliptic cone at angle of 
attack 

leeward surface pressure coefficient on elliptic cone at angle of 
attack 

differential operator 

interpolation constants in equation (Al) 
length of elliptic cone 
free- stream Mach number 


iii 



m slope of pressure distribution curve with respect to A* at 

horizontal plane of symmetry (sketch (b) ) 

n exponent characterizing pressure distribution at zero angle of 

attack 

n+ exponent characterizing windward incremental surface pressure 

distribution due to angle of attack 

n- exponent characterizing leeward incremental surface pressure 

distribution due to angle of attack 

S (n) normalized value of pressure distribution integral with respect to 

1 A* at angle of attack 

S^(n) normalized value of pressure distribution integral with respect to 
z/b at angle of attack 

S (n) normalized value of pressure distribution integral with respect to 

2 A* at zero angle of attack 

s(A*) normalized pressure distribution function (eq. (4)) 

y vertical coordinate 

z horizontal or spanwise coordinate 

a angle of attack 

AC r incremental chord-force coefficient due to angle of attack 

ACp incremental pressure on cone due to angle of attack in horizontal 

z plane of symmetry 

AC incremental surface pressure coefficient due to angle of attack 

(a+ windward surface, a- leeward surface) 

S c circular cone semiapex angle 

6 wo equivalent two-dimensional wedge angle, zero angle -of -attack case 

t 

6 wct equivalent two-dimensional wedge angle, angle -of -attack case 

Sy semiapex angle of elliptic cone in vertical plane of symmetry 

5 Z semiapex angle of elliptic cone in horizontal plane of symmetry 

6* proportionality constant (eq. (6)) 

y slope of correlation line (fig. 1) 

\p polar coordinate (sketch (a)) 
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A METHOD FOR PREDICTING PRESSURES ON ELLIPTIC CONES 


AT SUPERSONIC SPEEDS 
George E. Kaattari 
Ames Research Center 


SUMMARY 


A method is presented for estimating the pressure distribution over 
elliptic cones at supersonic Mach numbers and angle of attack. The method is 
based on an empirical correlation between experimental pressures in the sym- 
metry planes of elliptic cones and the pressures given by two-dimensional 
shock theory. The method is applicable for Mach numbers greater than 2, for 
cones whose ellipticity ratios range from 1 to 6, and for cones whose maximum 
semiapex angle is less than 30°. Results given by the method are shown to 
agree well with experimental values. 


INTRODUCTION 


Studies have shown that an elliptic cone can develop reasonably high 
lift-drag ratios and yet provide usable volume. Such a body, equipped with 
suitable stabilizing surfaces, offers possibilities as a supersonic flight 
vehicle. 

Exact theoretical solutions for pressure distributions on elliptic cones 
at angle of attack are difficult to derive and are usually numerical solutions 
to exact equations requiring machine computation. Approximate solutions for 
pressure distribution are given by Newtonian theory and by the tangent cone 
method. In many cases these approximate methods give pressure distributions 
that agree well with experimental values; however, there are ranges of Mach 
number and cone geometry in which good results are not obtained. Both methods 
fail to account for pressures over the leeward surface of a cone when it is 
shielded from the free stream by angle of attack. 

This report presents a semiempirical method with which pressure distribu- 
tions over the entire surface of an elliptic cone can be estimated with good 
accuracy for a wide range of Mach numbers, cone geometry, and angles of attack. 

ANALYSIS 


The analysis required for predicting pressure distributions on elliptic 
cones will be developed in two parts. First, the zero angle-of-attack case is 
considered. The analysis is based on an experimentally suggested invariance 
of chord- force coefficient with cone section ellipticity (ref. 1) and on a 


linear correlation of experimental pressures on the symmetry planes of ellip- 
tic cones with two-dimensional values from shock theory. This information is 
sufficient for predicting the pressure distribution over the cone surface. 
Next, the analysis will be extended to the angle-of-attack case. The pres- 
sures on the symmetry planes of the cone at angle of attack will also be shown 
to correlate in a simple way with values from two-dimensional shock theory. 
These symmetry plane pressures plus additional correlations to be described in 
detail then provide information sufficient to estimate pressure distributions 
closely on a cone surface at angle of attack. 


Elliptic Cones at Zero Angle of Attack 

Chord force - Experimental investigation of elliptic cones (ref. 1) indi- 
cates that for a given Mach number, cone length, and base area, the chord 
force is constant over a large range in the ellipticity ratio, b/a, of the 
cone cross section at zero angle of attack. This result is predicted by slen- 
der body theory in which the chord force depends only on the maximum cross- 
sectional area and is independent of shape. Newtonian theory also gives this 
result for small ellipticity ratio or, in any case, if the maximum angle of 
the cone apex is less than about 30°. The assumption was accordingly made 
that the zero angle-of-attack chord-force coefficient of an elliptic cone at 
a given Mach number is equal to that of a circular cone of the same length and 
base area. 

Correlation of conical and two-dimensional pressures - Typical correla- 
tions developed between experimental pressures (refs. 1,2) in the symmetry 
planes of elliptic cones and the corresponding two-dimensional or wedge pres- 
sures when the wedge angles are equal to the cone semiapex angles in the ver- 
tical, y, and horizontal, z, planes of symmetry are shown in figure 1. Linear 
correlations result at a given Mach number for a family of cones that have the 
same ratio of base area to length but differ in cross-section ellipticity. 

The circular cones are identified by filled symbols. The pressures at these 
coordinate points may be readily determined from supersonic handbooks (e.g., 
ref. 3). The appropriate cone for each family, fixed by the constant base 
area and length requirement, has the semiapex angle, 6 C , determined from 

tan 6 c = /(tan 6^) (tan 6~^) = /ab /l (1) 

The value of the slope, y = dCp/dCp^, wou ld then completely define the correla- 
tion line. The pressures in the symmetry planes of any elliptic cone in the 
family could then be determined from the correlation line with the two- 
dimensional (wedge) pressure coefficient corresponding to the cone angle in 
the plane of symmetry and the free-stream Mach number. 

Experimentally determined slopes y (refs. 4,5) including those of 
figure 1 and certain theoretical values (ref. 6) were plotted on logarithmic 
scales as a function of the ratio of circular cone pressure to wedge pressure. 
The general, linear correlation shown in figure 2 was found to result for the 
wide range of cone geometries and Mach numbers shown. This correlation is 
closely represented by 


2 



( 2 ) 


y = 



6 


The correlation lines of figure 1 are thus completely specified by C D /C^ 

F c p w 

and may be constructed for any family of elliptic cones without appeal to 
experiment since accurate theoretical values of pressure coefficients for 
circular cones and two-dimensional wedges are available. 



Sketch (a) 


Pressure distribution- A typical pres- 
sure distribution on an elliptic cone at 
zero angle of attack is represented in 
sketch (a) by the solid curve, Cp. The 
dashed curves represent a decompositon of the 

solid curve into Cr> and C n - curves. 

F zo Pzo P 

The abscissa. A*, represents the ratio of 
the area of the indicated sector to that of 
the whole quadrant. The variable A* is 
utilized since the integral of the pressure 
distribution curve, Cp, is then equal to the 
chord-force coefficient referred to the base 
area which, in view of previous discussion, 
is equal to the pressure coefficient, Cp c , 

of the equivalent circular cone. Symmetry 
considerations require that the derivatives 
of the pressure distribution curves with 
respect to A* be zero at A* = 0 and at 
A* = 1. 


The known pressures C^ , Cp , and C v are related as follows: 

Pyo Pzo 


>. -f 


C„ dA* = C T 


zo 


-f 


(C 


p ~ Cp) dA* - Cp - (Cp - Cp ) j s(A*) dA* 
^ zo r r zo r zo r yo 9 


>/■ 


or 


C - C 
Pzo P( 

Cp - C r 


Vo 4 


s(A*)dA* = S 2C (n) 


( 3 ) 

The term s(A*) represents the monotonical ly decreasing curve ( c p zQ - Cp) 

normalized to have the value of unity (1) at A* = 0. A simple, monotonically 
decreasing function having the required properties s(0) = 1, s(l) = 0, 
d[s(0)]/dA* = 0, and d[s(l)]/dA* = 0 was chosen to represent s(A*) 


S (A*) = 1 * c °^ ,(A * )n (4 

Values of the integral (eq. (3)) were found for various values of n. The 
results, represented by S^(n) as a function of n, are plotted in figure 3. 
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Since the pressure ratio 

( C P Z o " C Pc)/ (^zo " C Pyo) = S2c(n) 

is known, the value of n corresponding to S^fn) is also known from fig- 
ure 3. The pressure distribution curve, Cp, as a function of A* is then 
determined with equation (5) . 


c p - 


C n + C-n 
Pzo Pyo 


C P z o - C Pyo cos „(*.)» 


( 5 ) 


Elliptic Cones at Angle of Attack 

Pressures in the vertical plane of symmetry - The method for estimating 
the pressures on an elliptic cone in the vertical plane of symmetry at angle 
of attack is as follows: The required angle, 6 wa , for a two-dimensional wedge 

was determined so that the associated wedge pressure would equal the experi- 
mental pressure on the cone element in the vertical plane of symmetry. A 
linear relationship was assumed between the required angle, 6 wa , and the total 
inclination 6y + a of the cone element. The proportionality constant 6* 
was determined at zero angle of attack and assumed to be valid at angle of 
attack, that is 

6 6 
wo = = wa 

6y Sy + a 

or 

5 wa = 6*(6y + °0 



The pressure on the cone in the vertical plane of symmetry is then that of a 
two-dimensional wedge at the same Mach number and whose surface is inclined 
at the angle, 6 , as determined with equation (6). 

Equation (6) was verified by good correlation between predicted and 
experimental pressures as shown in figure 4 for a variety of cones at 
different Mach numbers and angles of attack. 


Pressures in the horizontal plane of symmetry - The method for estimating 
the pressure coefficient in the horizontal plane of symmetry on elliptic cones 
was based on an empirical modification of Newtonian theory. Newtonian theory 
gives the following expression for the variation of pressure in the horizontal 
plane of symmetry on an elliptic cone with angle of attack. 
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za 


C n cos 2 a 
Pzo 


or 




( 7 ) 


C Pzo C Pza _ AC Pz 
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Values for c p za given by exact theory in reference 7 for various circular 

cones, angles of attack, and Mach numbers were plotted in figure 5 with the 
coordinates indicated by equation (7). A linear correlation of the low angle 
(0° < a < 5°) data represented by solid symbols is evident. It was discovered 
that the ordinate, ACp^/sin 2 a, multiplied by the term, sin S z , brought 

into correlation the data of those cones whose "effective” slenderness param- 
eter, sin 6 , was less than unity. The higher angle-of-attack data approach 
the lower or Newtonian correlation. This is to be expected since at angles of 
attack approaching 90° the pressure, Cp , will be substantially zero by any 

theory, that is, ACp^ = Cp zQ . An interpolation technique for estimating the 

effect of angle of attack between the empirical low-angle correlation curve, 
AC pz /sin 2 a = 0.500 + 1.33Cp zQ , and the high-angle correlation curve, 

ACp^/sin 2 a = Cp^, is discussed in the appendix. The essential results of 

the interpolation are presented in figure 6 as A ^p z / c p zo as a function of 

angle of attack, a, for various values of the parameter, 

(0.500 + l-33Cp zo )/(Cp zo M w sin 6 Z ). 

Although the above correlation was based on the data for circular cones, 
it is assumed to be valid for elliptic cones as well. This assumption is 
indirectly confirmed by good chord-force predictions for highly elliptic cones 
at angle of attack. 



Sketch (b) 


Pressure distribution - Typical incremen- 
tal pressure distributions on an elliptic cone 
due to angle of attack are shown in sketch (b) 
Curves AC p ^ + anc * A ^p a represent the incre- 
mental pressure distributions on the windward 
and leeward surfaces of the cone, respectively 
In this case, symmetry considerations require 
that the slopes of these curves in the verti- 
cal plane of symmetry (A* = 0) are zero at all 
angles of attack. It is again assumed that 
these pressure distribution curves are of uni- 
formly decreasing or increasing value over the 
interval 0 < A* < 1. The curve of ^ C p ct+ has 

the slope -m at A* = 1 and the curve of ACp a _ 

has the slope +m at A* = 1. Simple monotonic 
functions of A* having zero slope at A* = 0 
and a nonzero slope at A* = 1 were used to 
represent the pressure distribution curves 
shown in sketch (b) . 


AC p a+ ' c p a+ ■ S? 


[(%>♦ ■ So) * - 4 S Z 

(8a) 
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AC p 


= C T 


- C r 


' '[(%» ' AC pJ cos [ ( ’ /2:itA * : ' n ] - 


ac t 


( 8 b) 


The average value of the integrals with respect to A* of equations ( 8 ) 
gives the incremental chord -force coefficient, &Cq, due to angle of attack as 
expressed in the following: 


AC C = 


[few - C Pyo) * *P,]W - [( C P yo - C P y „.) ; SC pJ S ic("-) 


yo fjya- 


- AC 


The functions S^fn) required in equation (9) are defined by the integral 


Pz 

( 9 ) 


■f 


S lC (n) = | cos[(Tr/2)(A*) n ]dA* 

A plot of S r (n) as a function of n is given in figure 7. 

1 Li 


( 10 ) 


The normal -force coefficient is the difference between pressure coeffi- 
cient curves ACp^ + and ACp^ (sketch (b)) integrated with respect to the nor- 
malized span variable, z/b = sin(iTA*/2). The result is 


C N = 


(c p - C p + AC p )s 1N (n+) + (C 
\ r ya+ r yo z / \ 


p ~ 

yo ya- z 


^ s i N (n- 


) 


( 11 ) 


it tan 6 


where the function S^fn) is defined by the integral 


s in w S J[ cos [(f) ( sin "‘ f) ]<• (f) (12) 

A plot of Sjjvj(n) as a function of n is given in figure 7. The normal-force 
coefficient is determined with equation ( 11 ) when the values for n are 
specified. 

The values for n required in equations ( 8 ), (9), and (11) are consid- 
ered next. Newtonian theory gives the result for the pressure derivative m 
indicated in sketch (b) that 


m 


dC 

— P- = (2 sin 2 
dA* V 


ir sin 2ct 
^z * 2 tan <5 

y 


(13) 
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fhe slope m of pressure distribution curves given by a more sophisti- 
cated (exact) theory of reference 7 for circular cones were compared with 
Newtonian values in figure 8. Surprisingly good agreement between Newtonian 
values and exact values was found as shown. 


The foregoing discussions in regard to incremental pressure distribution 
due to angle of attack and to the pressure derivative m in the horizontal 
plane of symmetry apply directly to circular cones that have a uniform pres- 
sure at zero angle of attack. In the case of elliptic cones, a modification 
to the Newtonian pressure derivative is required. This correction takes into 
account the fact that an elliptic cone does not have constant zero angle-of- 
attack reference pressure distribution and although the pressure derivative is 
zero at A* = 1, it rapidly assumes a value other than zero for A* < 1 . It 
was found for a wide range of cone ellipticity ratios that a good approxima- 
tion to the zero angle-of-attack pressure distribution was characterized by a 
curve whose average pressure derivative in the vicinity of A* = 1 is 
2(Cp - Cp ) as indicated in sketch (c) . Therefore, an improved slope for 

the incremental pressure distribution in the 
vicinity of A* = 1 is obviously given if the 
slope term 2(Cp zQ - Cp yo ) is added to the 

slope , m given by equation (13). The cor- 
rection to m should be subtracted for the 
leeward pressure distribution; however, a 
discontinuity in pressure slope would then be 
imposed at A* = 1. The corrective term 
2(Cp zo - c Pyo) was therefore applied in the 

positive sense to both windward and leeward 
pressure distribution calculations to avoid 
this discontinuity. The values for n+ and 
n- are related to this modified slope by 
differentiating equations (8) and recalling 
that dC p / dA* = 0 at A* = 1, 



Sketch (c) 


Curve AC n 

F a+ 




dC r 


a+ 


dA* 




^Pya+ ^Pyo 


+ AC P Z ) I (n+) 



or 



> C14a) 
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Curve AC p 

F a- 




or 


dCp 

r a 

cLA* 



m + 



(n-) 


2 [ m : 2 JSo - S o )] 

TT^Cp - Cp - ACp \ 

\ ^y o r ya- r z/ 



J 


(14b) 


Equations (14) with equation (13) give numerical values for n+ and n- . Pres- 
sure coefficient distributions, chord- and normal-force coefficients can then 
be calculated with equations (8), (9), and (11) and the curves of figure 7. 


COMPARISON BETWEEN EXPERIMENTAL AND PREDICTED RESULTS 


The comparisons between experimental and predicted results presented 
herein are not exhaustive but are typical of a larger body of such comparisons 
made during the course of this investigation. Direct comparisons between the 
predicted results of the present method and other simple methods were limited 
to the zero angle-of-attack case. 


Aerodynamic Coefficients 

The validity of the method for determining pressure distribution was 
assessed indirectly by comparison of the integrated pressure coefficients in 
the form of the aerodynamic coefficients, Cp and Cp with experimental values 
from reference 1 presented in figure 9. These experimental data represent a 
wide range in cone ellipticity, 1 < b/a < 6, at Mach numbers 2 and 3. The 
agreement between experimental and predicted values is considered good over the 
greater portion of the angle-of-attack range of the comparison. Discrepancies 
between experimental and predicted values are appreciable only in a few cases 
and then only when the angle of attack is greater than about 10°. This lack 
of good agreement at higher angles of attack is due to the effects of vortex 
flow in the wake. These effects are known to become significant with increas- 
ing angle of attack. The present method, based on inviscid flow, does not 
take into account any effects of vortex flow. 

It should be noted that a small experimentally determined skin-friction 
chord -force coefficient has been taken into account in calculating the pre- 
dicted values presented in figure 9 in order to make an equitable comparison. 
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Pressure Distributions 


Ze ro a ngle of attack - Experimental and predicted pressure distributions 
are presented for elliptic cones at zero angle of attack in figure 10. The 
predicted values include those of the tangent cone method and Newtonian theory 
as well as those of the present method. These pressure distributions indicate 
that the present method gives significantly better agreement with experimental 
values than do the other methods for cones of elliptic! ty ratio, b/a = 3 and 
higher. For lower ellipticity ratios, the pressure distributions given by the 
present method are as good or better than values given by the tangent cone 
method. The Newtonian theory values are consistently low. 

Ang le of attack - Experimental and predicted pressure distributions of two 
elliptic cones at angle of attack and at two Mach numbers are presented in 
figure 11. Predicted results of the other methods considered at zero angle of 
attack are omitted for sake of figure clarity. The agreement between experi- 
mental and predicted pressure distributions is considered satisfactory, partic- 
ularly at the higher Mach number of 6. A notable area of disagreement between 
predicted and experimental pressures exists- at the higher angle of attack 
(a = 20°) on the leeward surfaces of both cones at Mach number 3.09. This 
disagreement is due to the effects of vortex flow in the wake previously 
mentioned . 


CONCLUDING REMARKS 


An empirical method has been developed for predicting pressure distribu- 
tions over elliptic cones at angle of attack. The method is applicable for 
Mach numbers above 2, for cones whose ellipticity ratio, b/a, ranges from 1 to 
6, and for cones whose maximum semiapex angle does not exceed about 30°. 

The method is based, primarily, on an empirical correlation between 
experimental pressures on the symmetry planes of elliptic cones with two- 
dimensional values from shock theory; and, secondarily, on correlations of 
spanwise pressure derivatives with parameters suggested by Newtonian theory. 

Predicted pressure distributions and aerodynamic coefficients were com- 
pared with experimental values for air flows in the Mach number range from 2 
to 6, with cones of ellipticity ratio, b/a, ranging from 1 to 6, and angles of 
attack ranging from 0° to 20° . Satisfactory agreement between predicted and 
experimental values was found. The predictions given by the present method 
are generally in better agreement with experimental values than those given by 
the tangent cone method or by Newtonian theory. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, May 13, 1970 
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APPENDIX A 


DERIVATION OF THE INTERPOLATION FUNCTION USED IN PREPARING 
THE CURVES OF FIGURE 6 


As stated previously, it was found that the change in the pressure coef- 
ficient in the horizontal plane of symmetry at low angles of attack was 
expressed as ACp^/Cp^ = [(0-500 + i . 33Cp^^) / Cp^. q ] sin 2 a and at angles of 

attack near 90° as ACp^/Cp^ = s ^ n2 a - sketch (d) the slope of the 



Sketch (d) 


line 0-1 is (0.500 + 1 * 3 ^ C P Z o^^Pzo > re P" 

resenting the low-angle correlation curve. 
The slope of line 0-2 is 1, and represents 
the Newtonian and/or high-angle-of-attack 
correlation curve. The solid curve repre- 
sents a monotonic interpolation function 
of sin 2 a which is tangent to the low- 
and high-angle correlation curves in the 
neighborhoods of their respective validity 
(sin 2 a = 0 and 1.0). 


The simplest monotonic interpolation 
function was found to be a series expres- 
sion of three terms using the ordinate, 

AC d /Cr, , as the independent variable; 

* z ^zo 

that is. 


sin 


2 




(Al) 


In order to satisfy the known slopes and ordinates required of the interpola- 
tion function, the values of the coefficients, k, were found to be 


M sin 6 


ki = 


zo 


0.500 + 1.33C 


Pzo 


k 2 = -2(k x - 1), and k 3 = (ki - 1) 


(A2) 


Equation (Al) was evaluated for various numerical values for k for tabulat- 
ing sin 2 a as a function of (ACp^/Cp^) . This tabulation was used to con- 
struct curves of ACL, /CL, as a function of angle of attack, a, for various 

Fz Ezo 

values of the parameter, k^ . 
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Figure 4.- Experimental and estimated pressures in the vertical plane of 
symmetry of elliptic cones at angle of attack. 




Figure 5.- Effect of angle of attack on the pressure coefficient in the 
horizontal plane of symmetry of cones. 




Figure 7.- The functions S lC (n) and S lN (n). 
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Figure 8.- Correlation of exact theory (ref. 7) and Newtonian pressure 
coefficient derivatives, m, in the horizontal plane of symmetry of 
circular cones. 
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